Fluctuating initial conditions in hydrodynamics 
for two-particle correlations 

y. Hama^* R.P.G. Andrade^, F. Grassi^ and W.-L. QiaiP- 
^Institute dc Fisica, Universidade dc Sao Paulo, Brazil 

^Dcpartamcnto de Fisica, Universidade Federal de Ouro Preto, Ouro Preto-MG, Brazil 
DOI: will be assigned 



Event-by-event hydrodynamics, with fluctuating initial conditions, has shown to nicely 
reproduce several features of experimentally observed quantities in high-energy nuclear 
collisions. Here we discuss how it may help to understand, in a unified way, the various 
structures observed in the long-range two-particle correlations, both in nucleus-nucleus and 
p—p collisions. Suggestions of how experimentally this description could be tested are also 
discussed. 

1 Introduction 

In hydrodynamic approach of nuclear collisions, it is assumed that, after a complex process 
involving microscopic collisions of nuclear constituents, at a certain early instant a hot and 
dense matter is formed, which would be in local thermal equilibrium. This state is characterized 
by some initial conditions (IC) , usually parametrized as smooth distributions of thermodynamic 
quantities and four- velocity (see, for instance, [Hill). 

However, since our systems are small, important event-by-event fluctuations are expected in 
real collisions. Also, if the thermalization is verified at very early time, they should be very 
bumpy. In previous works, we introduced fluctuating IC in hydrodynamics [31 H], by using 
NEXUS event generator [5], and showed important effects on several observables. 

In this paper, we briefly survey some of the previous results [31 IHl [HI IH] and then discuss 
more recent results on long-range two-particle correlations. 

2 Some consequences of fluctuating initial conditions 

In Figure[l] we show the energy-density distribution in a typical event of fluctuating IC, gener- 
ated by NeXuS [S] for a central Au-fAu collision at 200A GcV. Observe that the distribution is 
very bumpy, as expected in real collisions, having a tubular structure in rj (leftover from initial 
particle collisions). 

Some consequences of such high-energy-density spots have been discussed in [6] . Because of 
high concentration of energy in small regions, each tube would suffer a violent explosion and, 
because of their small size, expand isotropically (in transverse directions). If such a tube is at 
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Figure 1: NEXUS Fluctuating Initial Conditions. 



the surface of the matter, certainly the outgoing part of this initial acceleration would remain, 
producing high-pT particles, which would be isotropically distributed in the momentum space. 
Thus, first we expect that high-p^ part of the px spectra is enhanced when fluctuating IC are 
used in our computations, in comparison with the results with averaged (smooth) IC. In the 
second place, we expect that the elliptic flow coefficient < V2 > suffers reduction as we go to 
high-pT region, due to the additional high-p^ isotropic components included now. A word of 
caution is necessary here. We are talking about effects of hot spots and not of the fluctuation, 
which makes V2 coefficient larger in more central windows because it makes the eccentricity 
bigger [7| . As for the rj dependence of V2 , we know that the average matter density decreases as 
\r]\ increases as reflected in the 77 distribution of charged particles, so when such a blob is formed 
in the large-jTyj regions, its effects appear more enhanced. Therefore, we expect considerable 
reduction of V2 in those regions. All these features have explicitly been verified in [6]. 

Another effect of small high-energy-density spots in the IC is manifested in the smaller HBT 
radii, as compared with the case of the smooth averaged IC. [9]. This has been shown both by 
using the Cooper-Frye prescription |IOj . and by the continuous-emission one 

Besides the effects of high-energy-density spots, fluctuations of IC imply evidently fluctua- 
tions of the resulting observable quantities. We have discussed that such fluctuations become 
quite large in the anisotropic- flow parameter V2 [HI [5], as has been effectively verified by 
experiments [T4] . 

3 Two-particle correlations in hydrodynamic approach 

One of the most striking results in relativistic heavy-ion collisions is the existence of structures 
in the two-particle correlations [El [TH [171 El [13 [2D] plotted as function of the pseudorapidity 
difference A77 and the angular spacing A^. The so-called ridge has a narrow Acj) located around 
zero and a long A77 extent. The other structure located opposite to the trigger has a single or 
double hump in A0; its Arj extent is not well established. More recently, the ridge structure 
has been observed also in pp collisions at LHC [21] . 

In an earlier work, |22| . we presented evidence that hydrodynamic approach reproduces all 
such structures in heavy- ion collisions. In |22| . the events computed by using the hydrodynamic 
code SPheRIO [4], starting from event- by- event fluctuating IC, generated by NeXus [5], were 
analyzed in a similar way to the experimental ones, in particular the ZYAM method was used 
to remove effects of elliptic flow. We later developed a different method to remove elliptic 
flow from our data and checked that all structures are indeed exhibited and other features well 
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reproduced (dependence on the trigger- or associated-particle transverse momentum, centrality, 
in-plane/out-of-planc trigger, etc) p3 l [24 l [25 l [26 ] . 

3.1 Mechanism of ridge formation - one-tube model 

As seen in Fig. [1] each NEXUS IC is very comphcated, so difficult to visualize how various 
structures in the two-particle correlations are generated. In order to clarify the origin of the 
ridge structures, we introduced in [23] a simplified model which would allow to follow closely the 
time development of the fluid in the vicinity of one of the high-cnergy-density tubes. Evidently, 
only those tubes located close to the surface of the hot matter can contribute to the correlations. 
Thus, in our simplified model, we replace the complex bulk of the hot matter by the average over 
many events, leaving just one typical tube close to the surface, like the one on the line 1 of Fig. [1] 
right. To simplify the computation, the longitudinal expansion is assumed boost-invariant and 
the transverse expansion is computed numerically (see details in [23]). 

Figure [2] shows the temporal evolution of the hot matter in this model. As seen, pressed by 
the violent expansion of the high-energy-density tube, the otherwise isotropic radial flow of the 
background is deflected and guided into two well defined directions, symmetrical with respect 
to the initial tube position. Notice that the flow is clearly non-radial in these regions. 




Figure 2: Temporal evolution of energy density for the one-tube model (times: 1.0, 3.5 and 8.5 
fm). Arrows indicate fluid velocity on the freeze out surface, thicker curve labeled by the freeze 
out temperature 0.14 GcV 

The resultant single-particle angular distributions for two different py-intervals are plotted 
in Fig. (3] left. As expected they show symmetrical two-peak structures. From this plot, we 
can easily guess how the two-particle angular correlation will be. The trigger particle is more 
likely to be in one of the two peaks. We first choose the left-hand side peak. The associated 
particle is more likely to be also in this peak i.e. with A(/) = or in the right-hand side peak 
with /S.(f) ~ +2. If we choose the trigger particle in the right-hand side peak, the associated 
particle is more likely to be also in this peak i.e. with Ac/) = or in the left-hand side peak 
with ^ — 2. So the final two particle angular correlation must have a large central peak 
at A0 = and two smaller peaks respectively at A0 ^ ±2. Figure [3] (right) shows that this 
is indeed the case. The peak at A(/) = corresponds to the near-side ridge and the peaks 
at A(/) ~ ±2 form the double-hump ridge. We have checked that this structure is robust by 
studying the effect of the height and shape of the background, initial velocity, height, radius 
and location of the tube [53] . 
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Figure 3: Angular distributions of particles in some different pT intervals (left) and resulting 
two-particle correlations (right), in the one-tube model. 

As stressed above, this simplified description, namely one-tube boost-invariant model, has 
been introduced just to clarify the mechanism of ridge-structure formation. However, it is 
remarkable that this simple model can describe well so many characteristics observed in experi- 
ments. For a more realistic simulation, we should consider more complex events such as NeXus 
events and average over the fluctuations. For an event like the one shown in Fig. [l] only the 
outer tubes need to be considered. The shape of the two-particle correlations for a single tube 
(in particular the peak spacing) is relatively independent of its features so the various tubes will 
contribute with rather similar two-peaks emission pattern at various angles in the single-particle 
angular distribution. For this single event, the two-particle correlation has a well-defined main 
structure similar to that of a single tube (Fig. [3]) surrounded by several other peaks and de- 
pressions due to trigger and associated particles coming from different tubes. These additional 
peaks and depressions have positions depending on the angle of the tubes between them. When 
averaging over many randomly fluctuating events these interference terms disappear and only 
the main one-tube like structure is left. The main advantage of this interpretation of ridge 
structures is that it involves essentially only the surface of the hot matter. The complexity of 
the kernel does not influence. 

3.2 Ridge in pp collisions 

Durig this Symposium, new results on two-particle correlations in pp collisions at LHC have 
been reported [27l [21] . Accordingly, high- multiplicity pp collisions show long-range {Arj > 4) 
nearside (A(/) ~ 0) ridge structure, which appears more enhanced when intermediate pt (1 < 
< 3 GeV) are considered. 

Such a structure can be produced also within our tube-based description. However, since 
NeXuS generator [5] is based on nuclcon-nucleon collisions as an ingredient, it can produce only 
one tube in the IC for pp collisions themselves. Evidently, this will not produce ridges. Some 
parton based generator is needed to a more complete description of pp data. 

Here, we shall show how such a structure can be produced within our description, by con- 
sidering a typical IC in pp, formed just by two tubes. Because of their small size, probably 
one-tube IC constitute the bulk of the fluctuating IC in high-energy pp collisions. However, 
although in a small fraction, there should be also events with two tubes produced in the initial 
collisions of constituents. In Fig. |4j left, we show such two-tube IC. K. Werner et al. discuss 
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nearside ridge, by using EPOS generator. [28], where the typical IC shown has a more eUiptical 
shape or two tubes with smaller distance. It is easy to understand that the two tubes begin 
expanding radially, collide with each other, driving the flow into two opposite directions, as 
shown in Fig. 21 middle, producing an azimuthal distribution of particles with two peaks sep- 
arated by TT. This flow is similar to the elliptic flow, but probably with sharper peaks. Final 
results of two-particle correlations in such events are two ridges, one at A(/) ~ and the other 
in the awayside at A0 ~ tt. In the CMS paper [5T], the awaysidc ridge is not mentioned. In 
our opinion, this occurred because it is really not visible due to the dominance of jet originated 
correlations. It would be nice to experimentally separate the ridge from the jet events, by 
measuring the three-particle correlations in the following way: 

Fix the first associated particle close to the trigger (A0i = (j>i ~ (j)t ^ Aiji = rji — rjt ^ 0). 
Presumably, we are picking jets in this case. On the other side, if we fix A0i ~ and A?7i > 
2, we are selecting ridge particles. Then, measure the correlation with a second associated 
particle in the (A(/)2 = 4>2 — 0t, A?72 = ir)2 — r]t) plane. We think that in the first case, we 
will see only typical jet correlations, with a high peak at Arj ^ A(p ~ and a high associated 
ridge at A02 ~ tt. And in the second case two low ridges separated by Acf) ^ tt, eventually 
with some contaminations of jet events. We call this kind of three-particle correlation 2+1 
correlation \24\ 126] . 




X[fm] X[fm] -3-2-10123 

Ai|> 



Figure 4: Time evolution of two-tube IC for pp at two different instants (r = 1 and 2.75 fm), 
and the resultant two-particle correlation as function of A(j). Here, the normalization has been 
done by using only one event, so the result should not be compared directly with data. 



4 Conclusions 

Hydrodynamic approach starting from event-by-event fluctuating initial conditions, with high- 
energy-density tubes, has shown to reproduce several features of experimentally observed quan- 
tities. 

With regard to the long-range two-particle correlations, it gives a unified picture for the 
nearside and awayside structures as observed experimentally. In heavy-ion collisions, a high- 
density tube located close to the surface of the hot matter divides the flow coming from inside 
into two currents, producing two-peak angular distribution. This two-peak distribution is the 
origin of both the nearside and the awayside ridges. In high-multiplicity pp collisions, probably 
there are a fraction of events with two high-energy tubes. Such IC produce azimuthal distribu- 
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tions with two peaks separated by tt, resulting in a two-particle correlation with two opposite 
ridges, one observed as the nearside ridge and the other hindered below the awayside structure 
associated with jets. 

See contribution by R. Andrade at this Symposium |29) for detailed discussion of the in- 
plane/out-of-plane effect within one-tube model. 
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